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Abstract 

Given integers ai, a2, . . . , a„, with ai + 02 + • • • + a„ > 1, a symmetrically constrained 
composition Ai + A2 + • • • + A„ = M of M into n nonnegative parts is one that satisfies 
each of the the n\ constraints { o-i^-Kii) ^ • tt G S'n }. We show how to compute 
the generating function of these compositions, combining methods from partition theory, 
permutation statistics, and lattice-point enumeration. 
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1 Introduction 



1.1 Constrained Compositions 



This work was inspired by the "constrained compositions" introduced by Andrews, Paule, 
and Riese in |^. We consider the problem of enumerating symmetrically constrained com- 
positions, that is, compositions of an integer M into n nonnegative parts 

M = Ai + A2 + --- + A„ = |A|, 

where the sequence (Ai, A2, . . . , A„) is constrained to satisfy a symmetric system of hnear 
inequahties. For example, the compositions Ai + A2 + A3 of M satisfying 

A7r(l) + ^7^(2) > A^(3) (1) 

for every permutation vr G S3, are known as integer-sided triangles of perimeter M 
114] . The number Am of incongruent triangles of perimeter M is given by 

J]AMg^= = (1-,2)(1_\3)(1_,4)- 

A/>0 Ai>A2>A3>0 ^ ^ ^ ^ ^ 

A2+A3>Ai 



However, 3 + 2 + 1 and 2 + 3 + 1 are different compositions (i.e., ordered partitions) of 
6 and counting the number A^^ of ordered solutions to ([1]) gives 

Ai+A3>A2 
A2+A3>Ai 

One could generalize this example in several ways. For example, moving to n dimensions, 
one could ask for the integer sequences (Ai, A2, . . . , A„) satisfying 

^tt{1) + '^7r(2) + • • • + A7r(„,_l) > ^T^ln) 

for all n! permutations vr € Sn- Another generalization would be to study, given positive 
integers k,l,m, the integer sequences (Ai, A2, . . . , A„) satisfying 

kK{i) + IK{2) > fnK{^) for all vr € Sn . 



Another example, considered in [3], was inspired by a Putnam exam problem [12^ Prob- 
lem B3]: Enumerate all compositions of M = Ai + A2 into two parts satisfying 

2Ai > A2 and 2A2 > Ai . 
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It is shown that 



2„,2 



Ai Aa l + xy + x^y 
y X ^ = — — , (3 

^ (1 - X?/2 (1 _ a;2y 



2A2>Ai 

giving a complete parametrization of all solutions. 

In [5], Andrews, Paule, and Riese demonstrate the suitability of the Omega package 
[3] for experimenting with problems of this sort and the power of Macmahon's partition 
analysis [2] to prove some elegant generalizations. 

The goal of this paper is (1) to formulate a generalization of the symmetrically con- 
strained compositions enumeration problem; (2) to show how this problem is connected 
to permutation statistics; (3) to show that the permutation statistics approach gives, for 
many cases, an effective computation method and, for certain cases, a way to derive com- 
pact formulas; and (4) to show that the insight provided by the geometry of lattice-point 
enumeration aids in the handling of the most general case. 



1.2 The Symmetrically Constrained Compositions Enumeration Problem 

Fix integers ai, 02, . . . , On- We are interested in enumerating compositions A = (Ai, A2, . . . , A„) € 
Z'"q that satisfy the n! homogeneous linear constraints 

aiK{i) + a2K{2) H \- (inKin) > for all vr G 5^ . 

Specifically, we are interested in computing the generating functions 

Fiz,,Z2,...,Zn) := J]z^4^•••z^ 
A 

and 

Fiq) := F{q, q,...,q)=Y, q^^+^-+-+>^- , 

A 

by exploiting the symmetry of the constraints. Note that because of the symmetry, there 
is no loss of generality in assuming that 

a-i < 0,2 < ■ ■ ■ < an , 

which we will do from now on. 

In Section 2, we show how to solve the enumeration problem when XlILi = 1. In 
certain special cases, we show that permutation statistics can be used to derive elegant 
formulas. We note that even this simple case is difficult for general purpose software like 
the Omega Package [3], Xin's improvement of Omega and LattE macchiato [10^ I13j. 
designed to enumerate solutions to linear Diophantine equations and inequalities. In Section 
3 we solve the general problem. We close this section with some notation and background 
on permutation statistics. 
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1.3 Permutation Statistics 



Throughout the paper, the following notation is used: [n]g = (1 — g")/(l — q)', [n]q\ = 
nr=i[^]g; and ia;q)n = U'i=oi'^-aq'). 

For a permutation vr = 7r(l)7r(2) • • • 7r(n) of [n] := {1, 2, . . . , n}, the descent set of it is 

D^ = {j: TT{j) >7r(j + l)}. 

The statistic des(7r) = j-D(7r)| is the number of descents of vr and the major index of vr is 
the sum of the descent positions: maj(7r) = ^- It is well known that 

(see, e.g., [H]). The joint distribution of des(7r) and maj(7r) over the set 5„ of all permuta- 
tions of [n] is given by Carlitz's g-Eulerian polynomial [6l[7|: 

n oo 

C^{x,q) = = Y[{1 - xq^)^[j]^ x^-\ 

Applying the definition of [j]q and the binomial theorem, we can rewrite this as 



Cn 



(1 — (/)" \i J 1 — q^x 



So, for example. 



Ci{x,q) = l 

C2{x,q) = l+xq (6) 
C3(x, q) = l + 2xq + 2xq^ + x^q^. 

If we take the limit as x — > g"" in ([5]) all terms except i = n in the sum are canceled by 
q)n+i = in the numerator, so 

Finally, for i < n — 1, let Sn^ be the set of permutations of {n\ that have no descent in 
positions {n — i, n — i + 1, n — 1}. Let 
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In [8], it is shown that 



^\xqn~k Cn-k{x,q) 



so, in particular, 



1 — xq 



2 Symmetrically Constrained Compositions when ^ = 1 
2.1 The Main Theorem 

Theorem 1. Given integers ai < 02 < • • • < fln satisfying '^^^i a-i = 1, i/ie generating 
function for those A G Z"q satisfying 

n 

ajX^^i^j^ > /or a// n £ Sn 



2S 



F(Z1,Z2,...,Z„)= ^ ^.^^^^^ J— 



where 



1 ifj = n, 

= ^ ifn>i>j>l, 

1 — (oi + • • • + Oj) if 1 < i < j < n. 

In particular, setting zi = ■ ■ ■ = = q yields 



(1- 9") n"=i'(i ^'-'^^ 

Proof. To simphfy notation, let 

F{z) = F{zi,Z2, ...,Zn). 

For 6 G Z", let 

Z — ^1 ^2 Z^ 

and for vr G 5„, let 

^TT = (^7r(l) ) ^7r(2) ) • • • ) ^7r(n) ) • 
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With 



7=1 



3 = 

we have 

F{z) = Y,z\ 
AeL 

Now we use the standard method of partitionmg the elements of L into classes L-j^ indexed 
by permutations it Sn- 

= I A G Z" : A^(i) > A^(2) > • • • > A7r(n), 

n 

aiX^(^i^ > for all a £ Sn, and 

1=1 

K{i) > K{i+i) if i G -Dtt I • 

Since the last condition guarantees that no A is in more than one class, L is the disjoint 
union 

Our goal now simplifies to computing 
because F (z) = X^^reSn ^'^ ^'^^ since 

K{1) > ^7r{2) > ■ ■ ■ > K[n) 

and since, by our assumption, ai < 02 < • • • < Om the n! constraints 

+ a2Acr(2) H ^ a„A^(„) > for ah G 5„ 

are all implied by the single constraint 

0-lK{l) + a2A^(2) H h anKin) > , 

so that we get the more compact description 

yn . -^7r(l) > K{2) > ■ ■ ■ > K(n) > and A^(j) > K{j+i) if J G -D^ 



i"^ ^ + 0-2K{2) H 1- 0.nK{n) > 

But this means that all look similar, except for the strict inequalities determined by D.,^. 
More precisely, if we let 

:= {A G Lid : Aj > Aj+i if j G D.„} 



and Gt^{z) := J2x^l then 

So it remains to find G-,^{z), the generating function for 

Ai > A2 > • • • > A„ > and \j > Xj+i if j € -Dtt 



for a given vr € 5„. 



aiAi + a2A2 + h OnAn 


> 




are given by the system 






-1 




ei 


1 -1 




62 




A > 




1 -1 




Cn-l 


«2 «3 • • • «n-l Q-ra 







(9) 



where 

^^fo ifj^D., 

We make use of the following lemma, a well known result in lattice-point enumeration. This 
version was formulated in [9] for easy application to partition and composition enumeration 
problems. 

Lemma 1. Let C = [cij] be an n x n matrix of integers such that = B = [bij] exists 
and bij are all nonnegative integers. Let ei,...,e„ be nonnegative integer constants. For 
each 1 < i < n, let Ci be the constraint 

Ci, lAl + Ci,2A2 H h Ci^nK > Gi . 

Let Sc be the set of nonnegative integer sequences A = (Ai, A2, . . . , A^) satisfying the con- 
straints Ci for all i, 1 < i < n. Then the generating function for Sc is: 



\&Sc 



j=lV-^l -^2 ) 



.A2 



Now let C be the matrix on the left side of ([9]). Then det(C) = ai + ■ ■ ■ + = 1, so C 
is invertible and B = has all integer entries: 

if j = n, 

hj = ( — (oi + • • • + Oj) a n > i > j > 1, 

if 1 < i < j < n. 
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If, in addition, ai + ■ ■ ■ + aj < for 1 < j < n — 1, the integer entries oi B = C are all 
nonnegative and Lemma 1 gives the generating function Gt^{z) and the theorem follows. 

To complete the proof, we show that if ai + • • • + a„ = 1 and ai < a2 < ■ ■ ■ < an, then 
for 1 < J < n — 1 we have ai + • • • + < 0. 

Let j be the smallest index satisfying 1 < j < n — 1 and ai + ■ ■ ■ + aj < 0, but 
ai + • • • + Oj+i > 0. Then Oj+i > —{ai + ■ ■ ■ + aj) > 0. Thus 

1 < aj+i < • • • < On , 

so 

1 = ai H h a„ > ai H h aj+i + n- j- l>n-j-l. 

So j = n — 1 and therefore ai -\- ■ ■ ■ + aj < for 1 < j < n — 1. □ 

In Section 2.3 we derive an algorithm based on Theorem 1 for efficient computation of 
F{q), given the Oj. In the next section, we give examples of how to combine Theorem 1 
with results on permutation statistics to derive formulas for F{q) in special cases. 



2.2 Applications 

Example 1 Given positive integers b and n > 2, let L be the set of nonnegative integer 
sequences A satisfying 

{nh-h+ 1)A^(„) > 6(A^(i) H h A^(„_i)) for aU vr € Sn • 

The case n = 2, 6 = 1 is the Putnam problem ([3]). Here a = [—6, —6, . . . , —6, n5 — 5 + 1], so 
by Theorem 1, 

p,. ^ n,.D. i'^''- ^ E.esM^'n-^'-''^ 
(1 - n"=i (1 - (1 - 'z") n"=i'(i - 9^+^"") ■ 

By dH), the numerator is just and simplifying gives 

1 — qn(nh+l) 



(1 - - q^b+i-^n 



This generating function was discovered by Andrews, Paule, and Riese and a complete 
parametrization was proved in [3] using partition analysis. 

Example 2 Given positive integers b and n > 2, let L be the set of nonnegative integer 
sequences A satisfying 

HK(2) + ■■■ + K{n-i)) >{nb-b- l)A^(i) for all tt e Sn ■ 
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The case n = 3, 6 = 1 is the integer-sided triangle problem ^ and the case n = 2, 6 = 2 is 
the Putnam problem ([3]). Here a = [—{nb — b — l),b,b, . . . ,b], so by Theorem 1, 

By 0, the numerator is 

f „n(bn—l) . „l—bn\ ( ^\(l—bn)n 
r l^n{bn-l) „l-bn\ _ W inj-gj' . 

^•^^^ 'g ~ ql^bn^n ' 

Simplifying further and dividing by the denominator gives 

1 _ fjn{nb-l) 

Fiq) 



(1 - - g^f*-!)" ' 

This generating function was also originally proved by Andrews, Paule, and Riese in [3]. 

Example 3 Given positive integers b and n > 2, let L be the set of nonnegative integer 
sequences A = (Ai, . . . , A„) satisfying the constraints 

{b + 1)A^(„) > 6A^(i) for ah vr G 5„ . 

The case n = 2, 6 = 1 is the Putnam problem ([3|). Here a = [—b, 0, 0, . . . , 0, 6 + 1] so by 
Theorem 1, 



By ([5]), the numerator is 

^bn. 



i=0 

Combining with the denominator and simplifying gives 



^ (l-g''")(l-g''"+") ^ (n\ _ {-q') 
(l-g")(l-g)" ^AiJl 



i=0 



bn+i 



Example 4 Given positive integers k < I, and n > 3, let m = + Z — 1 and let L be the 
set of nonnegative integer sequences A satisfying 

kK(n^i) + lK{n) > "zA^(i) for ah n e Sn ■ 

The case n = 3 and k = I = lis the integer-sided triangles Here a = [— m, 0, 0, . . . , 0, fc, /], 
so by Theorem 1, 

V TT n „j+mn-nk 
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Recall from ([8]) that Cn^ is the joint distribution of des and maj over all permutations with 
no descent in position n — 1. Then in F{q), we can split the sum over tt £ Sn into two sums 
according to whether or not i € D-,^. We get that the numerator can be written as: 

Using ([8]) and combining with the denominator gives (eventually) 

^ Cnjq''"', q){l - g"'-^) - Cn-l(g""', g)ng""'+»-l(l - g-"^) 

Details appear in [8]. 



2.3 Efficient Enumeration of Symmetrically Constrained Compositions 

We can compute the generating function F{q) for compositions satisfying the n\ constraints 

n 

«i A7r(i) > for all n £ Sn 



i=l 

via Theorem 1. The denominator is given explicitly, but the numerator is a sum of nl terms. 
However, regardless of the values of the Oj, the numerator of F{q), when simplified, is a 
polynomial with at most 2"""^ terms (one for each possible descent set). 

Let ui,U2, ■ ■ ■ be arbitrary and define polynomials G„ by 
We can compute G„ in the following way: Let 

where SSni is the set of all permutations it £ Sn that end with i. Then Gn = G^^^^ 

A permutation vr in SSni can be obtained uniquely from some permutation vf in Sn~i 
by replacing each j > i with j + 1 and then appending i at the end. The descent set of vr 
will be the same as the descent set of vf if the last entry of vf is less than i and the descent 
set of vr will be Djf U {n — 1} if the last entry of vf is greater than or equal to i. Thus we 
have the recurrence 

i—l 71—1 

with the initial condition G^-^^ = 1. We can simplify this a bit to get "Algorithm G": 

(i-i) 



G» = + (1 - Un-i)GllZl^ for i > 1 
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with d'^ = u„_iE"=iG'!^2i- 

Now, to compute the numerator of F{q) in Theorem 1 using Algorithm G, simply set 

Ui = g*-"("i-l ^"-i) for 1 < i < n and compute G^_^i\ 

If we use dynamic programming to implement the recurrence of Algorithm G, (e.g. 
"option remember" in Maple), then to compute Gn = G^^_^i \ at most O(n^) polynomials 
are computed. However, we must consider the time required to compute them. In order 

(i) 

to compute one of the , essentially we only need to add two polynomials. It is fair 
to assume that the time is proportional to the number of terms in the polynomials times 
the logarithm of the coefficient magnitude. So, overall, the time (and number of terms) 
grows roughly like 2" in the dimension n, but logarithmically in the coefficient size, which 
is considered polynomial time in fixed dimension. In practice, we found that we could 
compute F{q) for arbitrary a with ^ Oj = 1 within seconds for n < 11, in about 10 seconds 
for n = 12 and in less than a minute up to n = 15, using a naive implementation in Maple 
on a tablet PC running Windows XP. 

For comparison, there are existing software packages that, when given a collection of 
linear inequalities, produce the generating function for the integer points in the solution 
set. These packages include the Omega Package [3], Xin's speed-up of Omega [H], and 
LattE macchiato [10\ [13] . We used these programs to compute symmetrically constrained 
compositions in n dimensions, by giving as input the n! inequalities. The computation 
became infeasible when n > 4 for the Omega package and Xin's program. LattE was able 
to handle examples for n = 5 in under 10 seconds and n = 6 in under an hour. 

Thus exploiting the symmetry via Theorem 1 and Algorithm G makes a huge difference 
in what we can compute. 



3 The General Case 

3.1 A General Version of the Main Theorem 

We remove the requirement that ^ a-j = 1 and enumerate compositions A = (Ai, A2, . . . , A^) G 
Z"q that satisfy the n! constraints 

oiA^(i) + a2K(2) H ^ anK(n) > for all tt G 5„ (10) 

via the generating function F{z) = Ylx z^- 

Theorem 2. Given integers ai < a2 < ■ ■ ■ < an, with ai+a2 + - ■ ■+CLj < for 1 < j < n — 1 
and ai + a2 + ■ ■ ■ + an > 1, define the vectors Ai, A2, ■ ■ ■ , An G Z" as the columns of the 
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matrix 



and let 



Then 



02 H h On 


03 H h On 


04 + • • • + a„ 




On 


1 


-ai 


03 H h On 


04 + • • • + a„ 




an 


1 


-ai 


-oi — 02 


04 H + a„ 




an 


1 


-ai 


-oi — 02 


-oi -02-03 




an 


1 


-ai 


-ai — 02 


-oi - 02 - 03 




an 


1 


-ai 


-ai — 02 


-ai - 02 - 03 


■ ■ ■ -ai - 


■ ■ ■ — ttn-l 


1 




n 


f 














< Ci < 1 1 








i=i 


i=i 









F(z]= 'S^ >p iltgD^,pi=p,+ i 



where we take the product over all descent positions i o/vr for which the ith and the (i + l)st 
coordinate of p are the same. 

If, for some i, d divides every coordinate of Aj, we can replace Ai by Ai/d in Theorem [2] 
and thereby reduce the number of lattice points in P by a factor of d. 

Proof. The start of our proof is similar to that of Theorem [H except that we find it 
advantageous to view the compositions satisfying (flO]) as integer points in the cone 



K := i^x = {xi,X2, ...,Xn) 
From this point of view, 



n s 
G M>0 : X^ajX^(j) > for all vr G 5„ >. 
7=1 ^ 



The setup now continues in analogy with the proof of Theorem [H Like there, it suffices to 
study 

~ / 3:^1 > X2 > • • • > x„, > and xj > Xj+i if j G D.^ 

[ aiXi + 022:2 H h ttnXn > 

and the associated generating function GTr{z) := J2x(zk nZ'^ ^^^^ 
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First, we study the cone Kj^. The constraints of Ki^ are given by the system 



1 -1 
1 



ai a2 as • - 
and the inverse of the matrix on the left of is 



1 -1 

a?!— 1 a^j 



X > 0, 



(11) 



En 
1=1 



a2 + 1- a„ 


as + • 


• + a„ 


a4 + • 


• • + a„ • • 


an 


1 


-ai 


as + • 


■ + a,„ 


a4 + • 


• • + a„ 


an 


1 


-ai 


-ai 


- a2 


a4 + • 


• • + an 


an 


1 


-ai 


-ai 


- a2 


-ai - 


a2 - as 


an 


1 


-ai 


-ai 


- a2 


-ai - 


a2 - as 


an 


1 


-ai 


-ai 


- a2 


-ai - 


02 - as • • 


— ai — • • • — a„_i 


1 



The conditions on ai, 02, . . . , a„ guarantee that the inverse exists and that Ki^^ is a cone in 
M"q. Thus the columns Ai, A2, ■ ■ ■ , An of this matrix form a set of generators of Ki^ and 
by an easy tiling argument (see, e.g., [U Chapter 3]) 



i^id n z„ = I p + ^ CjAj -.per, c G z|o |; 



(12) 



in other words. 



E 



n-=i(i-4 

Before turning to K.,^, note that the generators Aj have a special form: we have 



A 



Ajj > Aj^ij for 1 < j < n and 
This follows because X]j=i aj > 1 implies that 

aj + Oj+i + 1- a„ > — ai — 02 - ■ ■ • — aj_i. 

Thus, if p € -fCjd n satisfies pj = Pj+i, then for any c G Z"q, if 

n 

r = {p + Aj) + ^c^Ai 



Ai+ij for j / i. 



2=1 



then 



rj > rj+i. 
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Now, what about i^^? It contains all points y G ETjd except those y with yi = yj+i for some 
i € -Dtt. By (fT2]) . every y € Kid has a unique representation as y = p + ^^j-^j fo'^ some 

c G Z>Q. Thus by the previous paragraph, yi = iff both pi = pj+i and Cj = 0. Now, in 
the same way as in p^ . 



n = < p + Cj^j : p € V, c G Z>0) and if j G and = Pj+i then Cj > > 

j=l J&Dt,, pj=pj+i 



A,: 



pePnZ" nj=i ~ 



Thus 



□ 

In the special case of Theorem 1, X]"=i o-j = 1 and the origin is the only lattice point in V. 
3.2 Efficient Computation for the General Case 

Give ai < ■ ■ ■ < an with Oj > 1, once we find the generators Ai, A2, ■ ■ ■ , An, and 

the lattice points in V, we can again use Algorithm G to efficiently compute F{q): By 
Theorem 2, setting z = {q,q, . . . ,q), 



where = xi + ■ ■ ■ + Xn for an ?i-dimensional vector x. 

The denominator is easy. To find the numerator, for each point p G V Ci Z", set 



Ui 



q\M if p. = p._^_^ 
1 otherwise 



and then compute G^*^^^ ■ 



Now the running time also depends on |PnZ'"|. This can grow linearly with the magni- 
tude of the entries (rather than the logarithm of the magnitude), even in fixed dimension. 
However, when iPriZ"! is of moderate size, this computation method can be quite effective. 
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